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ABSTRACT 


‘A digital computer program for the geometrically 
nonlinear analysis of totally arbitrarily loaded 
shells of revolution (SATANS-II) was modified to more 
accurately account for the conditions at the pole, of 
the shell. This program was used to determine the 
buckling load of shallow Spherical sheils of various 
sizes when subjected to static axisymmetric, dynamic 
axisymmetric, and dynamic nearly axisymmetric 
step- pressure loads of infinite duration. A 
comparison was made between the new buckling results 
and previous results obtained without the new pole 
routine. The comparison, revealed a significant change.’ 
in the buckling pressures, due solely to the change in 
the pole routine. The new static axisymmetric, 
dynamic axisymmetric, and even the dynamic asymmetric 
critical buckling pressure loads appear to be fairly 
reliable results for perfect, shallow shells. 
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the modulus of elasticity of the shell 

the rise of the spherical cap at the pole 
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the mass density of the shell 
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the nondimensional critical pressure 
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sphere 
a column matrix containing the coefficients 


of the a term in the series expansion of the 
applied load 
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the normal distance from the axis to the cap. 
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directions, respectively 
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of the shell 

the nondimensional time 

the time 


a reference time 
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the displacements in the s, © and 
directions, respectively 

nondimensional series coefficients of U, V, W 
a nondimensional measure of the volume of 

the shell deformation 

the peak in the time history of the 


parameter ¥ 
the displacement in the ¥ direction 


: th , 
in the n harmonic 
the nondimensional time increment 


distance between stations 


the nondimensional parameter governing the 
magaitude of the load applied in the asymmetric 
harmonics . 

the coordinate aormal to the surface of the shell 
the circumferential angle measured about 

the axis of revolution 

a nondimensional geometric parameter used 

to describe the spherical cap 
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the normal distance from the base plane to 


the middle surface of the undeformed cap 


ACKNOWLEDGEMENT 


I wish to acknowledge the help that I have received in 
the course of my thesis research. 

I am particularly indebted to Professor Robert E. Ball, 
who provided me with the guidance and aid necessary to 
complete my project. 

I am also grateful to Professor Johann Arbocz of 
CALTECH, who provided the initial deck of cards from which I 
worked. 

I also wish to give credit to the NPS programming 
consultants, and in particular Richard Donat, for the expert 
assistance he provided in debugging my changes’ to the 
computer program. 

I especially want to express my sincere appreciation to 
my wife, Roberta, for her patience and assistance in typing 


of this thesis. 


In 1973 a digital computer study was presented by Ball 
and Burt {1] for the dynamic buckling load of clamped 
shallow spherical shells subjected to axisymmetric and 
nearly axisymmetric step-pressure loads. A static buckling 
analysis of the same spherical shells had been carried out 
in 1970 by Stilwell and Ball [2]. In these two studies the 
digital computer program SATANS~I [3] was used to calculate 
the critical buckling pressures for a large range of shell 
sizes. Other studies of the puckling of shallow shells have 
been conducted by Huang {4,5], by Stephens and Fulton [6], 
by Lock et al. [7], by Stricklin [8], and most rcecentiy by 
Akkas {9]. In Reference 1 the results from these other 


studies, except for those by Akkas, are compared with the 


results from SATANS-I for both static and dynamic buckling. 
*In the axisymmetric static analysis tae comparison with the 
results obtained by Huang [4] revealed that the SATANS-I 
results were higher than Huang's results for several sheli 
sizes. In the dynamic, axisymmetric buckling analysis the 
SATANS-I results again either agreed closely with, or were 
somewhat higher than , the results by Huang [5], Stephens, . 
and Fulton (6]J, and Stricklin [8]. However, it was noted 
then that there was a general lack of consistent agreement 
among any of the sets of results. AS a consequence, it 
appeared at that time that the axisymmetric buckling problem 
had not yet been totally resolved and that additional 
studies would be appropriate. 


In the asymmetric dynamic buckling analysis of Reference 
1 the few comparisons that could be made for the critical 
load also indicatec that the SATANS-I results may be too 


high. A ccmparison of the recent estimates for the 
asymmetric dynamic buckling load obtained by Akkas [9] with 
the SATANS-I results also reveals the SATANS-I results to be 
well above those of Akkas [9]. However, it should be noted 
that the results obtained by Akkas were from his attempt to 
obtain a lower bound on the critical asymmetric load. This 
bound on the buckling load is obtained without the execution 
of a complete transient response analysis on the asymmetric 
Part of the response of the shell, as is done in SATANS-I. 
In Akkas*' analysis (Problem 1) the transient nonlinear 
axisymmetric response is computed, and a determinant is 
examined for possible bifurcation into asymmetric noticn at 
each time step. The minimum load at which the determinant 
becomes zero is defined as the lower bounce of the. critical 
load. 


As a consequence of the generally high buckling loads 
predicted by SATANS-I, a re-examination of the static and 
dynamic buckling of the shallow spherical shell was made in 
an attempt to determine the possible cause, or causes, of 
the high buckling loads. In our search we discovered that a 
nodification of the manner in which the pole conditions are 
numerically approximated significantly lowered the buckling 
loads to values that are now in good agreement with the 
other results. The new procedure for handiing the pole 
condition is given in section III of this thesis. The new 
buckling results are given in section V. 


In addition to the pole condition modifications and the 
new buckling results the author has also made another 
Significant change to the SATANS family of codes. In 
particular, the SATANS-II program for the geometrically 
nonlinear analysis of totally arbitrarily loaded shells of 
revolution, developed by Ryan [10] in 1972 to handle more 
complex and larger problems, was modified to make the 
computer memory requirement a variable quantity. This 
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quantity is specified by the user to fit the particular 
problem being run. It eliminates the large core requirement 
of SATANS-II for small problems and allows for much larger 
problems to be solved than could be solved by SATANS-II. 
The new program with the pole condition and semory 
modifications will hereafter be called SATANS-IIA. It is 
a@escribed in section II. 


SATANS-II was developed by Ryan [10] from SATANS-I and 
incorporated the full trigonometric expansion of the applied 
load and solution vector, and introduced the handling of 
imperfections into the code. These modifications allow the 
analysis of shells under totally arbitrary loads, as well as 
imperfection studies on actual shells with measured 
imperfections [11]. Unfortunately, the original deck of 
cards for SATANS-II was destroyed. Professor Johann Arbocz 
of CALTECH had a listing of SATANS-II and punched a deck of 
cards with the changes to SATANS-I given in that listing. A 
copy of this deck was sent to Professor Ball. These cards 
have been added by the author to the original SATANS-I 
described by Ryan [10] and a complete version of SATANS-II 
has been reconstructed. SATANS-IIA is a modification by the 
author of the reconstructed SATANS-II program. A listing of 
SATANS-IIA can be found in Appendix A. The listing contains 
an example problem for the dynamic analysis of a clamped, 
truncated cone subjected to an impulsive loading which is 
uniform along the meridian and varies in a cosine 
distribution over one-half of the circumference. This 
problem is a sample problem “suggested by the Lockeed 
Missiles and Space Corp. [12]. A condensed version of the 
output from the example problem is given in Appendix B. 
Input data preparation for SATANS-IIA can be found in 
Appendix C. The basic users manual, which includes 
preparation of input subroutines and the theory of the 
program, is contained in Reference 3, which can be obtained 
through COSMIC (M70-10098, LAR-10736), or ASIAC [13]. A 
users manual which includes preparation and handling of 
imperfection data within the SATANS programs can be found in 
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Ref. {10]. The above information, along with the followiag 
discussion, will inform the user on the capabilities and 
proper use of SATANS-IIA. 


The modification of the SATANS-II program to make its 
core requirement variable was accomplished by putting in a 
single dimension statement at the beginning of the progran, 
with subsequent dimensioning within the subroutines to only 
the first element of the vector or matrix. This isa 
convenient feature of the FORTRAN-IV language in which the 
program is written. The actual vector and matrix sizes are 
transmitted to the subroutines by an individual parameter 
list. Construction of the initial dimension statement and 
core request size is as follows: 


The basic size of the program on the IBM-360/67 Digital 
Computer, without the initial dimension statement, is 
272,000 bytes. This figure includes approximately 19,000 
bytes of buffer space required for execution. Within the 
Main dimension statement are fifteen variables. However, 
only three parameters are needed to specify the sizes of 
these fifteen variables. 


Let a= The number of stations along the meridian of the 
shell times the number of harmonics considered. 
Let b= a, plus two fictitious stations times the number 


of harmonics considered. 
Let c= The number of harmonics considered. 


The main dimension statement would then be constructed 
as, 
DIMENSION P(4,4,a), DEE(4,4,a) - DST(4,4,a),- X(4,a), 
PHIXB(a), PHITB(a), Z(4,b), Z0(4,b), 
Z2(4,b), 23(4,b), ZDOT(4,b), IS(99,c), 
JS(99,c), ID(99,c), JD(99,c) 
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The 99's above limit 
run and an unlinited 
core requirement for the 


the user to 99 karmonics in any one 
number of meridional stations. The 
general case would be, 


272,000 + 216a + 80b + 1584c = bytes of core required. 


For a sample calculation of the core requirements 
consider the example of a spherical cap with 40 stations 
along the meridian, and an asymmetric analysis with two 
harmonics. Therefore, 


a= 40(stations) x 2(harmonics) = 80 
b= 80 + 2 x 2(harmonics) = 84 
c= 2 (harmonics) 


Thus, for the variables P, DEE, DST, 
3 x (4x4x80) = 3840 (words) x 4 = 15,360 bytes 


for the variable X, 


4 x (80) = 320 (words) x 4 = 1280 bytes 
for the varibles PHIXB, PHITB, 
2 x (80) = 160 (words) x 4 = 640 bytes 


for the varibles Z, ZO, Z2, 23, ZDOT, 
5 x (4x84) = 1680 (words) x 4 = 6720 bytes 


lastly, for the varibles ID, JD, IS, JS, 
4 x (99x2) = 792 (words) x 4 = 3168 bytes 


Therefor, the total size of the main dimension 
statement would be 27,168 bytes. This figure would be 
rounded up to the nearest even thousand bytes, ie. 28,000 
bytes. Finally, the core requirement for this example 
problem would be 


272,900 + 28,000 = 300,000 bytes. 


| 
| 
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The SATANS code is based upon Sander's geometrically 
nonlinear equations under the conditions of small strains 
and moderately smail rotations. The formulation is in four 
second order nonlinear partial differential equations in 
terms of U, V, WwW, and ae where U, V, and W are the 


meridional, circumferential and normal displacements 

respectively, and #M is the meridional bending moment. The 
s 

neplinear partial differential equations in the coordinates 


Ss, @ , and t are reduced to uncoupled sets of linear 
differential equations in s and t by expanding the variables 
in trigonometric series in the circumferential coordinate 9, 
and treating the nonlinear terms as pseudo loads. The first 
and second derivatives in the meridional coordinate s are 
replaced by the conventional central finite difference 
approximations, ie. 


en 
N 

~. 
i] 


1/2M ( {Zz}. ~ {Zp ot (1) 
itt i-1 


and 


2 
{zy = I/O ({z},.. - 2 {2}, + {2}, .) (2) 
1 1+ 1 pt 


1 1 


where {z}. is the vector of U, V, W, and M at the i 
he s 


station, 4 is the uniform dimension between stations, and 


primes denote partial derivatives with respect to s. 
Applying these approximations to the governing set of domain 
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equations leads to 


fe), 4... * ES] &), + (4), 2)... = 153, (3) 
i i-1 i i i i+ 1 


When the shell does not have a pole, fictitious stations 
one increment off of the shell are introduced at each end. 
Both the governing domain equations and the boundary 
conditions are applied at the two boundary points. Thus, 
all finite difference approximations to the derivatives, 
including those of the boundary conditions, are of order 


2 
2 However, prior to the development of SATANS-iia, the 


treatment of the conditions to be applied at a poie at 
either end of a shell was handled by a simple Euier forward 
or backward difference approximation to the first 
derivative, with truncation error of order 4. For example, 
for a pole at s= 0, where i= 1, the first derivative at the 
pole was approximated with 


(zF, 21/8 (iB * E} ye (4) 


At the time this procedure for handling the pole 
conditions was developed (1967) it was thought that this 
would not significantly alter the solution. However, it has 
since been discovered that such is not the case. 


For the new pole routine, an expanded forward difference 


approximation of order A : is used at s= 0 which takes into 
account the two stations after the pole, instead of just one 
Station after the pole as in the Euler schene. This 
approximation is 


Y 21/24 (3 4 - ‘ 5 
{2}, / ( =, * {2}, {2} .) (5) 


16 


The conditions to be imposed upon the dependeat 
variables at a pole are derived in Reference 14. They are: 


For N= 0, u =v =w’ =m’ = 0. 
4 i 2 st 


Applying equation (5), these conditions can be put into the 


matrix forn 


- 


1000 u 0000 u 0000 u 
0100 v + 0000 v + 0000 v = 0. 
0 0-3 0 w 004 0 wW 0 0-1 0 “ 
00 0-3 m 000 4 D 00 0-1 n 

s’ 1 s‘ 2 sj 3 


where the above 3 matrices are DL, DG, and DF within the 


SATANS programs. 
Por N= 1, u tv =uw=w=nmn =0, 
where the plus sign applies at an initial pole, and the 


Minus sign at a final pole. The matrix fcrm for these 


conditions is 


-3 000 u 4000 u -1 000 1 
141 0 0 v + 0000 v + 0000 Vv = 0. 
0010 w 0000 wW 0000 Ww 
000 1 a 000 0 m 0000 n 
s' 1 s s* 3 
For N=2, u= v= w= m = 0 
s 
the matrix form is 
1000 u 0000 u 0000 
0100 v + 0000 Vv + 0000 v = 0. 
0010 w 0000 w 00q0q0 4 
00 0-3 n 000 4 m 00 0-1 m 
BS 4 Ss” 2 Ss 3 
For N > 2, u= v= w= n= 0 
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and DL= identity matrix, DG= DF= null matrices. 


The sclution procedure in SATANS is an elimination 
scheme and starts with 


z} =- (P z + {x 6 
Gy. cep i, {x}, ’ (6) 
where the values in CP), based upon the Euler approximation 
are defined in Reference 14. The higher order approximation 
defines a hew bey: This new ane is obtained by 
Simultaneously solving the pole conditions 
[DL] {2}, + [DG] {2}, + (DF] {2}, = {0}, (7) 
and the domain equation at station 2 next to the pole 
Cc Zz + (B z + [A z = 8 
[Cc], (23, Nee ae (Alt), {g} ’ (8) 
to eliminate a Thus, 
sy = fay" - [c {z} - [B] {z . 9 
GI. cal, ({93 ced i, Bd 35) (9) 
Substituting equation (9) into equation (7) gives 
4 
DL Z + {DG Zz + [DP A =aliG Zz = 
[ Per, [ ] 3, [ PLAY, ( ta), Led th, 
B z = 0. 19) 
C5), i.) ( 
Combining like coefficients of the {z} vector leads to 
( (DL) - [DF] [AJ*{[c])) {z} + DG] - [DF] [Aa}'EB 
{DL} -[ } fart 1, {2}, ( (DG) - (DF) [ i 1) 


5 = = a e 1 
(23, (DF ] tay. {3}, (11) 


Finally, solving for (2) | yields 
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- 


Z2}. =- [DL - DF x at xc DG - DF x A’ x B z 
{ }, [ A 2! ( 2 3? { 1, 


js . - he . 12 
+ { DL - DF x a x o1 [ - DF x ve) {9}, (12) 


~1 
Thus P = - [DL - DF x A* xc DG - DF x A’ xB and 
al }, f > ci] a 2 3! 


-1 
x}. = (DL = DF x a' xc - DF x a" - The new 
{ ,, [ 9 od C 5! {3}, 
{ nee matriy has been placed into the "PMATRX" subroutine of 
SATANS-IIA and the new iy, vector has been placed in the 


"FORCE" subroutine. 


A listing of the pole routine may be found in Appendix 
D. To inccrporate this new routine into a SATANS-I or-II 
program, first proceed to the "PMATRX" subroutine and remove 
the fifteen cards that are between, but not including, 
"IF (NN.GT.2) GO TO 90" and "11 CONTINUE". These cards are 
located after statement number "14" and just before 
Statement number "11". Replace the cards removed by the 
ones listed in Appendix D which read from" C IN PMATRX " to 
" 90 M3=MN". Then proceed to the "FORCE" subroutine and 
remove statement number "10". Replace statement number "10" 
with the nine cards listed in Appendix D which read from" Cc 
IN FORCE" to "DO 11 I= 1,4" +. Also place " COMHON 
/IBL5/IBCINL, IBCFNL " into the common area of the " FORCE " 
subroutine. 


This completes the implementation of the new pole 
routine into either SATANS-I or II. 


IV. PROBLEM DESCRIPTION 


The geometry of the shallow spherical shell used in this 
study is identical to that used in Reference 1. Briefly, 
the shallow shell can be specified by the non~dimensional 


parameter X , where 


2 #174 1/2 
A=2 3 (1-V )j (ah) (1) 


H is the rise of the shell, h is the thickness, and Y is 
Poisson's ratio. The mass density of the shell is ao. All 
shells analyzed had the following dimensions; 


Radii of Curvature R = Be = 250 inches 
s 

Thickness h = 0.25 inches 

Modulus of Elasticity E = 30,000,000 psi 

Poisson's Ratio y=0.3 


All buckling pressures obtained will be listed as a 
percent of the classical buckling pressure of a complete 


sphere, q4° where 


2 2 172 
a= C2E ( we Tk al Ge es aa le ea (2) 


Forty stations were used over the meridian. The 
nondimensional time increment § t, where 


2 1/2 
ee ry % — / &) ’ (3) 


was taken as 0.05 for 3000 time steps, which is a total 
nondimensional time of 150. In addition, the axisymmetric 
analysis was repeated with a larger time step of 5 t= 0.2 
for a total time of 600. In this study m was selected such 
that t is equal to T. The necessity for the long response 
time is explained in Reference 6. 


In the axisymmetric analysis only the N= 0 harmonic is 
considered. However, in the asymmetric analysis a_ second 
harmonic is excited by applying an incremental load in that 
harmonic. In addition, analyses of the shells X = 6, 7.5, 
and 1% were made using five harmonics. The step pressure- 
load for the axisymmetric harmonic is 


I id q) {1}, (4) 


and the step pressure load for the asymmetric second 
harmonic is 


(n) (n) 
Pe Ee, oe Peis (5) 


{¢ 


(n) 
where n> 0, and € is taken as 0.0001. The value taken 


for the second harmonic in the asymmetric analysis was the 
Same aS the critical harmonic for the static buckling 
analysis presented by Stilwell and Ball [2]. When there was 
an uncertainty as to which was the critical static harmonic 
the two harmonics in question were both tested. Run times 
using SATANS-IIA with a two-harmonic analysis for 3000 time 
steps and 40 stations on the meridian took an average of 28 
Minutes on the IBM 360/67. 


The parameter used to determine the minimum load at 
which dynamic buckling occurs is the peak value of ¥, called 
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q » where ¥ is defined as 


MAX 


«= ¢*_ . (0 A 
v= cow dr / e = ar (6) 
C) ° 
ris the normal distance from the axis to the shell, a is 


(0) 


the maximum value of r, W is the normal displacement of 


the axisymmetric response and $€ is the vertical distance 


from the kase plane to the undeformed shell. The ¥is a 
measure of the volume of the shell deformation. The Fortran 


statements computing V and ae are given in Appendix E. 


When working a problem that requires these calculations the 


nineteen cards are inserted directly into the "DYNAMIC" 
subroutine right after the "IF" statement that calls tke 
"OUTPUT" subroutine. 


For convenience, the response in each asymmetric 


0) 
harmonic is also measured using equation (6), with ra 


: (n) ® 
replaced with W - The parameter V for the asymmetric 


harmonics does not represent a volume of deformation as it 


does for the axisymmetric harmonic. It can, however, be 
used to indicate the relative excitation of the asymmetric 


harmonics. 


The buckling criterion for both the axisymmetric and the 
asymmetric dynamic buckling analysis defines the critical 
load as that load P where a very small increase in P causes 


‘a very large increase in - : This is the same criterion 


as that used in Ref. {1]. 
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Ve. RESULTS AND DISCUSSION 
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A. STATIC AXISYMMETRIC BUCKLING ANALYSIS 


Table I presents the new results from the _ static 
axisymmetric buckling analyses for >» = 4 through 13 using 
the new pole routine. The two upper curves in Figure 1 
present a comparison of the new results obtained by 


SATANS-IIA with those obtained by Stilwell and Ball [2] 


using the SATANS-I program. As can be seen in this figure, 
fairly significant changes in the buckling load occurred in 
the neighborhood of » = 4,5, and 9; and somewhat smaller 
differences cccurred in the region A = 10 through 13. The 
upper data points in Figure 2 present the comparison of the 
new results from SATANS-IIA with those obtained by Huang 
{4}. This comparison shows a very good agreement between 
the two sets of results, except for the largest values of Xd. 
The new results have eliminated the differences that existed 
between the SATANS-I results and Huang's results. 


B. DYNAMIC AXISYMMETRIC BUCKLING ANALYSIS 


Figure 3 presents the new results for the peak value of 


ve versus P for the various values of » tested. Table 


It presents all of the new results for the dynamic 


axisymmetric buckling load. These loads are selected from 


figures constructed just like Figure 3. In every case, 


except for » = 4, a value of P slightly above the a 


value caused a ee indicative of buckling, as well as a 


nonconvergence of the iterative solution procedure. 


The lower two curves of Figure 1 present a comparison 
versus > of the new axisymmetric dynamic buckling results 
with the previous buckling results obtained by Ball and Burt 
{1]. In every case the new critical pressure is lower than 
the critical pressure obtained using the Euler approximation 
at the pole. 


The lower data points of Figure 2 present a comparison 
of the new results with those obtained by Huang [5], by 
Stephens and Fulton [6], and by Stricklin [8]. Just as in 
the case of the static axisymmetric buckling analysis, the 
new results compare much more favorably with the other 
results than did the results of Reference 1. Its 
interesting to note that the new results now tend to be 
slightly lower than the other results, whereas the results 
of Reference 1 were aigher for almost all values of AX. 


C. DYNAMIC ASYMMETRIC BUCKLING ANALYSIS 


Table III presents the new results for the critical 
pressures obtained from the dynamic asymmetric analysis. 
The second harmonics, or critical static harmonics, used in 
the analyses are also presented in Table III. A comparison 
of the critical pressures from the asymmetric analyses, 
Table III, with the critical pressures from the axisymmetric 
analyses, Table II, reveals that only the shell A =6 
buckled at a load below the axisymmetric buckling load. For 
the shell 2A = 7 the critical buckling load was slightiy 
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larger when asymmetric motion was considered. In all other 
cases the buckling was not influenced by the presence of the 
second harmonic. These new buckling results and those by 
Ball and Burt [2] are plotted in Figure 4 The new results 
can be seen to be significantly different from the SATANS-I 
results, where the asymmetric buckling loads were lower than 
the axisymmetric loads for five out of the ten values of 
tested. 


Except for » = 6 and 7, the relationship between ae 


and P for the N= 0 harmonic, in the two-harmonic analyses, 


was found to be essentially identical to the relationship 
found in the axisymmetric buckling analysis shown in Figure 
36 Table IV A presents the re versus P data for both the 


N= 0 harmonic and the second harmonic, for all values of A 


tested, -except for X = 6. Note that, except for X=7 and 


11, the for the asymmetric harmonic is generally very 


small, even when the ae for the N= 0 harmonic indicates 
that the shell has buckled. Thus, except for the shells A= 


6 and 7, the presence of the asymmetric motion does not 
influence the axisymmetric motion, and except for the shells 
» = 6, 7 and 11 the asymmetric motion is very small prior 
to buckling in the axisymmetric harmonic. 


A more detailed analysis of the shell dX = 6 has been 
conducted since it was the only shell that revealed any 
Significant axisymmetric sensitivity to asymmetric motion. 
This shell yas studied using two two-harmonic anaiyses (N= 
0, 1 and N= 0, 2) and a five-harmonic analysis (N= 0, 1, 2, 
3, and 4). Figure 5 and Tables IV B and IV C contain values 
of ie versus P for both of the asymmetric harmonics, N= 1 
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and N= 2, in the two two-harmonic analyses, as well as the 


values of ae for the axisymmetric harmonic, N= 0. Figure 


6 and Table IV D present the values of ene versus P for the 


N= 0,1,2,3, and 4 harmonics from the five-harmonic study. A 
comparison of the critical buckling load predicted from the 
results of the two two-harmonic analyses in Figure 5 with 
the critical load from the five-harmonic analysis obtained 
from Figure 6 shows that the presence of the additional 
harmonics results in the shell buckling at a slightly lower 
load (0.50), with significant motion in the N= 1 harmonic 
instead of the N= 2 harmonic (see the norconverged solution 
at P= 0.51), which is the critical harmonic for static 
asymmetric buckling. Studies using five harmonics have also 
been conducted for A = 7.5 and X= 11. As can be seen in 
Table IV D the critical harmonic for X= 7.5 remained N= 3; 
however, significant motion occurred in that harmonic at P= 
-41 and .44. In the case of A = 11, relatively large 
asymmetric sotion occurred in the asymmetric mode of N= 5 
vice 6 at a value of P= .46. 


The comparison of the new results for the critical 
pressure for dynamic asymmetric buckling with those obtained 
analytically by Stricklin (8], by Akkas [9], and 
experimentally by Lock et al ({7] is iliustrated in Figure 7. 
The comparison reveals an agreement with Stricklin in every 
case, in general a higher value of Y nage than those obtained 


by Akkas, and most importantly a very good agreement with 
Lock's experimental results. 


When making the comparison between the new results and 
those obtained by Akkas, it is necessary to look at the 
differences in the problem solution parameters used in the 
two studies. For example, buckling results obtained froma 


SATANS-IIA using the same time increment as used by Akkas, 


§ t= .2 for 3000 time steps, were significantly higher 
than those using the time step of §t= .05 for many values 
of A. Furthermore, the new results had , in some cases, 
instances of buckling occurring as far out in time as 130. 
Akkas, to shorten computer run times, observed the cap only 
for a time of less than 5. Furthermore, only the harmonics 
N= 1 or 2 or 3 were studied by Akkas for shells A= 5 
through 12. If the critical harmonic is not studied, the 
predicted load will be too high. Thus, it appears’ that 
Akkas' lower bound loads may not be true lower bounds. 


Two additional features of the shell response should be 
noted. First, shells X» = 6, 7.5, and 11 exhibited a 
non-buckled response in the axisymmetric harmonic to a load 
larger than the defined critical buckling load. fhis can be 
seen in Tables IV A and IV C. Second, and most importantly, 
the buckling load proposed by Ball and Burt [1], and used 
here, defines. buckling to occur when the Le in the 


axisymmetric harmonic undergoes a large change due to a 


small change in P. Another criterion for dynamic buckling 
in the asymmetric analysis discussed in Reference 1 is to 
define the buckling load as that threshold load that 
initiates significant growth in the asymmetric harmonic. 
Re~examination of the ae versus P data in Table IV A 


through D reveals that shells A= 6, 7, and 11 exhibited 


relatively large asymmetric motion at loads smaller that the 


defined buckling load when compared with other te values 
AX 


for those shells, even though the numbers themselves were 
\ 


small when compared with the axisymmetric harmonic. ShellsaA 
27.5 and 12 appear to be borderline cases. If the alternate 
criterion for buckling is used, the critical buckling loads 
for shells >= 6, 7, and 11 become 0.47, 0.45, and 0.45, 
respectively. The shells » = 7.5 and 12 could have buckling 
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loads as low as 0.40 and 0.44, respectively. These values 
are more conservative than the definition based upon 
axisymmetric response. These five shells are the same five 
shells that exhibited an asymmetric buckling load lower than 
the axisymmetric buckling load in Reference 1. 
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Vi. UMMARY AND CONCLUSIONS 


a Ss == Sle Se eS 


A digital computer program for the geometrically 
nonlinear. analysis of totally arbitrarily loaded shells of 
revolution (SATANS-II) was modified to more accurately 
account for the conditions at the pole of the shell. This 
program, called SATANS-IIA, was used to determine the 
buckling load of shallow spherical shells of various sizes 
when subjected to static axisymmetric, dynamic axisymmetric, 
and dynamic nearly axisymmetric step-pressure loads of 
infinite duration. The cap sizes ranged from A= 4 to 13 
including \ = 7.5. A comparison was made between the new 
buckling results with the improved pole handling routine and 
the results that did not have the new pole routine. The 
comparison revealed a significant change in buckling 
pressures, due solely to the change from an order A _ finite 
difference approximation of the first derivatives at the 


pole to an approximation of order & . These new critical 


pressures ate in very good agreement with the results from 


other studies of the same spherical shells. This good 
agreement with other results, which came about as a resuit 
of the modification of the pole handling routine, is a 
strong indication that the manner in which the pole 
condition is handled is vital to the accuracy of the 
solutions obtained. 


In the asymmetric analysis, two harmonics were included 
for most of the shells; the axisymmetric harmonic and one 
asymmetric harmonic. Five-harmonic analyses were conducted 
for three of the Shells. Two buckling criteria for the 
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asymmetric aralysis were considered. One defined buckling 
as that threshold load that caused a large increase ina 


deformation farameter, Rs in the axisymmetric harmonic. 


The other, more conservative than the first, defined 
buckling as that threshold load that caused a large increase 


in the ee value for the asymmetric harmonic. Both values 


have been presented. 


The new static axisymmetric, dynamic axisymmetric, and 
even the dynamic asymmetric critical buckling pressure loads 
appear to te fairly reliable results for perfect, shallow 
shells. The effect of realistic imperfections remains to be 
determined. 
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Figure 3 - PEAK DEFLECTION VERSUS P, AXISYMMETRIC AND 
ASYMMETRIC CASES FOR VARIOUS VALUES OF A (SATANS~ITA) 
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Figure 5 - PEAK DEFLECTION VERSUS P FOR THE ASYMMETRIC 
ANALYSES OF 4 = 6 (N=0,1 AND N=0,2) 
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Figure 6 - PEAK’ DEFLECTION VERSUS P FOR THE ASYMMETRIC 


ANALYSES OF A = 6 (N=0,1,2,3,AND4, ONLY N=0,1,AND2 PLOTTED) 
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E IV Dynamic asymme 


1. TABLE IV A. Two-harmonic analyses for all values of 
dX exept A= 6. 
A = 5 N= 0 and 2 N= 0 and 1 
ee Oe ee | ae ee ee 
ju= 0 }.1659 | .1676 _|_.6606 {| w= 0| 1675 | .6606 | 
z_2 J.aoos7e7 [.aooanse | oper as a feona348 | 0897, 
ae Sy ea 


[= o[ .7653 
[a= 1 [0010382 _] 
X = 7, N= 0 and 3 


a Ce ee ee ee ee 

is o|ssese [sven overa[1oos [wah ast foe 

~04323 
53 

fe i 

w= 3] 05997 | .06252 | .03809] 


} = 7.5, N= 0 and 3 


CE OS SP RA 
= o[-0703 | -07290[ -o7va9 [ .2636[ area [2074 
ae 


N= 0] .200 
N= 3] 0003276 | 


db = 8, N= 0 and 4 
w= 0| .05893 | .0607 |.0s20 | 1960 | i709] 4957 
a eee! 

mmo] 2197 
naa] 0000326 


w= o[ .03239[ 03347 
IN=_5}.0000281]. 0000472 |.00004125|.00002103|.0000449|.000114 


\} = 11, N= 0 and 6 
Ce rT PT a aT | 50] 
F.004595 


ha=_6 | 
A = 12, N= 0 and 7 


IN= 0{  .03236 . 03316 . 08633 
-00004214 | .0004561 | .00005158 


A = 13, N= 0 and 8 


| 
| .02119 | £02185 -06637 | .07844 07381 
- 00001148 |.00001134 |. 000006607|.000008245|.000119 | 
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2. TABLE IV B. Two-harmonic analyses with N= 0 and 1, 


A = 6 cnly. 
[ep | eve | 49 | 50 515253 
[w= o| .1350 | .1385 | .1421_ | 1460 | . 1499 | . 4453 | 
[y= 1].0002797|.000195 | .000245 | . 000926 


TABLE IV C. Two-harmonic analyses with N= 0 and 2, 


Tren] es | ee] a? | ee | | 50 
000239 


3. 
A = 6 only. 


ie 
ran 
I el oe ee 


TABLE IV D Five-harmonic analyses for selected 


4. 
snells. 


X = 6 N= 0,1,2,3, and 4 
Fae ea Ma a 

1347 |.1382_ | .1460[.3797_| 1498 | . 2200 | 

N= +].00021 |.02108 |.02215 | .003676|.3743 | .01276 | .0616 | 
N= 2|.0187 |.0953_ |.1069 |.1507 |.1502_ |.1279 |. 2385 | 
N= 3|.000181|.006237|.01437 |.00163 |.o405 |.0123 |. 03978] 
N= 4|.0031 _|.04757 |.05428 |.ou402 |.o5a96 |.ou9s |. 064 
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A = 7.5 N= 0,1,2,3,4 and 4 
eC ee ee ee ee 


-00004855; .00004198 |.00006093 |.0002737 |.0001167 
-0001164 00007456 |.0004184 |{.0000982 }.000788 |.0003188 


N 


= 
u 


cS (Sea a 
[03920 [.04008 [0099 [ouss | coast | - 1000 
[-0005755].001263 | 0009774 L.00136a | 002657 | .001565 
009568 | 01508 
erry 
302068. 


.009568 |.0140 | .0124 02239 02759 
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APPENDIX A 


LISTING OF SATANS-IIA 
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